Quantum Hall Effect of Massless Dirac Fermions in a Vanishing Magnetic Field 
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The effect of strong long-range disorder on tire quantization of tiie Hall conductivity Oxy in 
graphene is studied numerically. It is shown that increasing Landau-level mixing progressively 
destroys all plateaus in a^y except the plateaus at a^y ~ =Fe^/2/i (per valley and per spin). The 
critical state at the Dirac point is robust to strong disorder and belongs to the universality class 
of the conventional plateau transitions in the integer quantum Hall effect. We propose that the 
breaking of time-reversal symmetry by ripples in graphene can realize this quantum critical point 
in a vanishing magnetic field. 

PACS numbers: 72.10.-d,73.21.-b,73.50.Fq 
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Graphene, an isolated layer of graphite [H, displays 
a remarkable quantization of the Hall conductivity 
when subjected to a magnetic field. As shown in Refs. [3 
and , a^y is measured to be either a negative or a pos- 
itive half-integer n 1/2 in units of Ae'^/h. This should 
be contrasted with the usual integer quantum Hall ef- 
fect (IQHE) for which quantum Hall plateaus of the Hall 
conductance in SiMOS-FET or in GaAs/AlGaAs hetcro- 
junctions occur at positive integer values n in units of 



In the vicinity of the two non-equivalent cor- 



ners K and K' of the first Brillouin zone of graphene, 
the non- interacting electronic dispersion is linear, i.e., 
at the so-called Dirac point it realizes a noninteracting 
massless Dirac Hamiltonian. As the Hall conductivity 
of a massive two-component Dirac fermion in (2 + 1)- 
dimensional space and time has long been known to be 
a half- integer in units oie^/h 0,0], the observed quan- 
tization of the Hall conductivity of graphene can be at- 
tributed to 4 independent flavors of two-component Dirac 
fermions with the spin-degeneracy accounting for 2 and 
the valley-degeneracy accounting for another 2 flavors. 

This explanation fails to account for Anderson localiza- 
tion due to disorder in graphene, as recently emphasized 
m Ref. 0. Disorder can introduce both inter-valley and 
intra-valley scattering between the Bloch states in the 
valleys centered at K and K' 0, [l^l ■ The physics of 
localization for the Bloch (Landau) states of graphene 
is predicted to depend sensitively on the nature of the 
disorder and, in particular, on the range of the spatial 
correlations of the impurities. When disorder is short- 
rangcd; inter-valley scattering has large matrix elements. 
In this case, all Bloch states are localized at zero magnetic 
field @], while the IQHE is wiped out by a sufficiently 
strong disorder 0. On the other hand, when disorder is 
sufficiently long-ranged, inter-valley scattering is negligi- 
ble, and the metallic phase remains stable to disorder if 



time-reversal symmetry is present [11|, |12|, [13|, [IJ] . 

As for the conventional IQHE, the quantization of the 
Hall conductivity in graphene requires the existence of at 
least one critical state in each impurity-broaden Landau 
level (LL) induced by a magnetic field. The destruction 
of the quantum Hall plateaus in the conventional IQHE 
with increasing disorder strength [isl . [l6j can be under- 
stood in terms of a critical state in a disorder-broaden LL 
migrating (levitating) to higher filling fraction as soon 
as the disorder induces significant LL mixing [17]. The 
same phenomenon has been shown to be operative for all 
impurity-broaden LLs in graphene when the disorder is 
short-ranged [§]. 

In this letter, we argue on the basis of numerical calcu- 
lation that, when the disorder is sufficiently long-ranged, 
all but one critical states in graphene undergo levitation 
as the disorder strength increases. The exception is the 
state at the Dirac point that remains critical whatever 
the disorder strength is. Consequently, graphene with 
long-range disorder has the remarkable property that the 
quantum Hall plateaus with a^y = ±e^/2/i per spin and 
per valley survive in the limit of strong disorder. In the 
language of the renormalization group (RG), the scaling 
flows of the longitudinal {uxx) and transverse {uxy) con- 
ductivities for graphene with strong infinite-range dis- 
order are determined by two attractive fixed points at 
(T^y ~ zLe'^ /2h and one repulsive fixed point at a^y = 0. 
For comparison, the corresponding phase diagram for 
graphene with strong short-range disorder is conventional 
in that it is characterized by a single attractive fixed 
point, the insulating phase with a^^ = i^xy — [sl- fioj. 
Of course, the range of the disorder is always finite in a 
sample of graphene [l3| . However, the characteristic scat- 
tering length induced by short-range impurities could be 
longer than the phase coherence length for some range 
of temperature. Indeed it has been found, on a sample 
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of graphene with a high-mobility (/i ~ 2 x lO^cm^/Vs), 
that the longitudinal conductivity at the Dirac point at 
zero magnetic field never falls with decreasing temper- 
ature down to 10 mK [l^; i.e., intervalley scattering is 
irrelevant in this range of temperature. 

The honeycomb lattice of graphene has two atoms per 
unit cell on sites labeled A and B. Linearization of the 
noninteracting lattice Hamiltonian in the clean limit, ig- 
noring the spin degrees of freedom, yields [H 

n (^^ , Hk = vpo- ■ [-ihV + eA] (1) 

where Tix' is the transpose of Ti-K, Vp is the Fermi ve- 
locity at the Dirac point, and the Pauli matrices a act 
on the sublatticc degrees of freedom. In a finite mag- 
netic field B = V X A = (0, 0, B) the spectrum of Hk 
consists of LLs with eigenvalues £"„ — sgn{n)fvujQyJ\ri\ 
with n Gil [2^. The energy and length scales are set by 
uoq = \/2vpIIb and Ib = \/h/eB. The filling fraction 
;^ > (i' < 0) is the ratio of the number of occupied 
(empty) states above (below) the Dirac point to the LL 
degeneracy. 

In this work, as we neglect inter-valley and spin- 
dependent scattering channels, every LL is four-fold de- 
generate. The problem of Anderson localization that 
we shall consider is defined by the random Dirac equa- 
tion \Hk + yif )] V' = £^/' with the static vector poten- 
tial in the Landau gauge A = (0, Ex, 0) and V a ran- 
dom static scalar potential. We choose two-dimensional 
space to be a square with area L^. We introduce a 
cutoff Tic = 9 for the LLs. The disorder potential 
is V(r) = Ej^r^tjcxp(-|r - Rj\'^ /2(f) / {2i:(f) with 
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FIG. 1: (color online) Disorder-average Hall conductivity a^y 
as a function of the filling fraction v for the dimensionless 
disorder strength (a) V /hujo = 0.4, (b) 0.7, and (c) 1.1. Upper 
inset: Disorder-average density of states (DOS). Lower inset: 
Size dependence of a^y at T/hujo = 1.1. 



r G R^, d a length scale, and u/d? = \uj\/(P an energy 
scale. The signs in Uj = ±u and the scattering centers 
Rj are chosen randomly and independently at A^imp lo- 
cations. The disorder strength is estimated to be the 
broadening of the LLs FV4 = 2TTu^Nin-,p/[ie% + 2(P)L^] 
within the self-consistent Born approximation [l^. In 
the following we set d = Q.TIb- 

We evaluate the transverse conductivity with the Kubo 
formula 

he^ f li-n[{m\v^\n)(n\vy\m)] \ 

"^^^"^..<t<..' ^' 

where \n) denotes an eigenstate with its eigenvalue e„ of 
Hamiltonian Hk + the corresponding velocity op- 

erator is V := vpcr, and the over line represents disorder 
averaging. The dependence of the disorder-average Hall 
conductivity a^y on the filling fraction is shown in Fig. [T] 
for three values of the dimensionless disorder strength. 
For small LL mixing induced by the disorder, wc observe 
that the curve (a) has plateaus centered at half-integer 
values and with rounded corners up to v = ±(nc 4- 1/2). 
There are fewer well-defined Hall plateaus at moderately 
strong disorder strength as shown in the curve (b). Only 
the plateaus centered at v = ±1/2 are seen in the curve 
(c) when the disorder is strong relative to hujQ. The size 
dependence of <7xy at strong disorder (c) shown in the 
lower inset of Fig. [1] is consistent with the expectation 
that axy has a step at = in the L — > 00 limit. The 
quantization of a^y / {e'^ / h) for all half-integers larger in 
magnitude than 1/2 correlates with an oscillatory de- 
pendence on the energy of the disorder-average density 
of states (DOS) as is illustrated with the upper inset of 
Fig. [TJ This correlation is not present for the plateaus at 
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FIG. 2: The Thouless conductance gr as a function of the 
filling fraction u when (a) T/hcuo = 0.4, (b) 0.7, and (c) 1.1, 
for various system sizes, L/£b = 15 (□), 25 (o), and 35 (V). 
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V = ±1/2. LL mixing induced by the disorder is thus re- 
sponsible for the disappearance of all Hall plateaus with 
Wxy \ > while the plateaus aX v ^ ±1/2 are robust 

to strong disorder. 

We now turn our attention to the dissipative com- 
ponent of electrical transport by studying the Thouless 
number qt = (|Ae|)/fc. The Thouless energy Ae is the 
energy shift induced on energy eigenvalues by using an- 
tiperiodic instead of periodic boundary conditions [l^. 
The mean level spacing is 5e = l/(L^p(e)) with p[e) 
the DOS. The angular brackets denote the typical value 
(Ae) = cxp(ln Ae). The dependence of gr on v at differ- 
ent values of L/£b is shown in Fig. [2] for weak (a), mod- 
erate (b), and strong (c) disorder strength. We observe a 
decrease of the Thouless number with the deviations of 
the filling fraction from the integer values n, which is well 
correlated with the oscillatory dependence of the DOS in 
Fig. [T] This decrease gets more pronounced the closer n 
is to n = 0. Local minima of the Thouless number are 
reached at half-integer filling fractions. At fixed filling 
fraction, generally decreases as the system size L/^b 
increases, while only at the center of an impurity-broaden 
LL it is independent of L/£b. This is consistent with the 
existence of a critical state close to integer filling frac- 
tions. The decrease of gets more pronounced when 
approaching half-integer filling fractions. Increasing dis- 
order strength as shown in Figs. [2ljb) and[2l^c) weakens 
the dependence of qt on both v and Lj Ib around all half- 
integer filling fractions larger than 1/2. A decrease of qt 
remains clearly visible at v = \/2 for strong disorder. 

For any filling fraction, a localization length ^ is ex- 
tracted from the fit griL) = gQCxp{—L/^) and 1/^ is 
plotted in Fig. for weak, moderate, and strong dis- 
order. Figure[3]Ja) is consistent with an oscillatory depen- 
dence of ^ on J/ that correlates with the DOS in Fig. [T] 
and also shows that ^ increases for v > 1 with disorder 
strength r/(ftci;o). The latter feature is familar in the 
conventional IQHE, where the LL mixing makes the lo- 
calization weaker in the high-field regime [l^. In each 
impurity-broadened LL of Fig. El^a), there exists a v for 
which ^ takes a maximum and much larger value than 
L. This u corresponds to a local maximum of g^ and 
to a critical single-particle state. In the limit of no LL 
mixing, critical states are found at integer values of v. 
Figure Elb) shows that all but one critical states are 
increasing functions of \v\ as the amount of LL mix- 
ing induced by the disorder increases. The exception 
is the v — Q critical state that remains pinned to the 
Dirac point for all disorder strength in Fig. [3l[b). Con- 
sequently in the weak magnetic field limit fnoQ/T — > 
all single-particle states in the hole (electron) regime are 
localized: (Jxy/ie^/h) = —1/2 (+1/2) with the step dis- 
continuity in axy occurring at i/ = 0. This behavior is 
different from the one for the conventional IQHE systems 
where ax y ~ » (Hall insulator) in the weak magnetic field 
limit 1^, 1^ . It is also different from the one observed in 
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FIG. 3: (color online) (a) Dependence on v of the inverse 
localization length 1/^. (b) The phase diagram for the quan- 
tum Hall effect of massless Dirac fermions. For each impurity 
broaden LL, the phase boundaries are identified from the L- 
independent qt (circles). Half-integers indicate —(Jxy/{e^ /h). 
(c) Widths of the half-maxima of the Thouless conductance 
(closed symbols) and the slope of the Hall conductivity (open 
symbols) as functions of L. (d) Scaling flow in cTxy — gr plane: 
Oxy and Qt at L/Ib ~ 15 (□), 25 (o), and 35 (V) are plotted. 



numerical studies of the effects of short-range correlated 
disorder on the IQHE in graphene which reported 
strong localization behavior in the vicinity of = over 
a wide range of magnetic field. This difference can be 
used to distinguish whether the dominant source of dis- 
order in graphene is short- or long-range. For v close 
to I' = and from weak to strong disorder. Fig. [3fc) 
is consistent with the scaling Ansatz axy = fxy{L^^°'"v) 
and gT = fxxiL^^""i^) with Ui, « 2.3. Correspondingly, 
1/C ~ W\"'% l/(ATi^) ^ L"" where I^tv is the half- 
width of the peak of at = 0, and daxy/dv ^ L'^" 
with a„ = I /ay. Moreover, gr = fxx{fxy{crxy)) must 
hold close to = 0, a prediction confirmed by the scaling 
flow in Fig. [3Jd), where {axy^gr) computed at L = 15, 
25, and 35 are seen to move towards the attractive fixed 
points at (±e^/2/i, 0) with increasing L. 

The phase boundary separating the quantum Hall 
phases axy = —e^/2h from axy = ±e^/2/i in Fig. [Hb) 
extends to the zero magnetic field limit fiwo/T —>■ 0. This 
strongly suggests that the Dirac point is never localized 
by long-range disorder, however strong the disorder is. 
The Dirac point of graphene with long-range disorder re- 
mains the critical point of the IQHE plateau transition 
for any weak magnetic field. 

We close this letter by studying the effects of effec- 
tive time-reversal symmetry breaking random perturba- 
tion on the zero-magnetic field limit Hloq/T — s- of the 
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FIG. 4: (color online) The /^-functions for the random Dirac 
Hamiltonians Q; l3{g) is independent of the mean value of 
V. 



Dirac Hamiltonian ([T]). The random singlc-particlc Dirac 
Hamiltonian is [6] 

n :=cT ■ [-ifiwFV + a(r)] + 'm{r)az + V{r), (3) 



where a.(r), m{r), 
tributed with the 
,5^^A^exp(-q2g;2/2)J(q 

(y,ai,a2,TO) whereby 



and V{r) are 
covariances 



Gaussian dis- 



Av 



fM)fM') 

if (/^=0,l,2,3) 

for ^ = 0, A^ = 



for /.t = 1,2, and A^ = Am for ^ = 3. The length 
scale d is the range of the potential that varies up 
to 1/30 of the minimal system size. The "random 
vector potential" a(r) describes the effect of ripples in 



corrugated graphene [2l|, 122] ■ Ludwig et al 



in Ref. S 

conjectured that the random Dirac Hamiltonian ^ 
flows into the critical point describing the quantum Hall 
plateau transition. Using the numerical method from 
Ref. [l3|, we present in Fig. |3] the beta function (3{g) 
for the Kubo conductivity g = (Jxx/ {^^ /h) which we fit 
against I3{g) = -l/(27r^5^) + £> exp(— 27r(;) Q whenever 
time-reversal symmetry is broken {—iayl-C^iay ^ Ti.) by 
either A^ > or Aa/ > 0. The fitting parameter D is 
optimized when D ^ 7. Figure [4] is consistent with a 
fixed point at g ~ 0.6 which is attractive in g. The beta 
function in the limit Aa,m with fixed Ay > is 
always positive (open circles in Fig. U]), as reported in 
Refs. [13[ and [1J|. This limit is equivalent to the limit 
hujo/T ^ in Eq. ([1]). Both limits induce a crossover 
from the unitary class with a topological term to the 



symplectic class with a topological term [ll|, |12| . This is 
different from the limit Hloq/T — > in the conventional 
IQHE that realizes a crossover from the unitary class 
with a topological term to the orthogonal class. The 
fact that the orthogonal class is always insulating in two 
dimensions thus causes all critical states to levitate to 
large ly as Hujo/T 



We conclude by pointing out the interesting possibility 
that ripples in graphene alone can realize scaling flows to 
the IQHE plateau transition point even in the absence of 
a magnetic field, i.e., graphene with long-range disorder 
is very close to the quantum critical point of the IQHE. 
This scenario is consistent with recent low-temperature 
transport measurements of high-mobility graphene sam- 
ples that imply a temperature independent minimal con- 
ductivity 
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